Proton superconductivity and the masses of neutron stars 
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The unexpected temperature evolution of the neutron star in the Cassiopeia A supernova remnant 
(Cas A, for short) has renewed tremendous interest in the cooling mechanisms of neutron stars. In 
particular, the formation of superconducting protons and superfluid neutrons deep inside the cores 
of neutron stars have become focal points of the discussion. The purpose of this letter is to add a 
new aspect to this discussion, which focuses on the connection between proton superconductivity 
and the masses of neutron stars. 

Assuming (as is currently the case) that the temperature evolution of Cas A is largely controlled by 
superconducting protons, we study a series of phenomenological proton-pairing models to determine 
how deep into the stellar core superconducting protons actually penetrate. This allows us to establish 
a heretofore unknown relationship between the mass of the neutron star in Cas A and the penetration 
depth of the superconducting proton phase. This relationship can be used to either predict the depth 
of the superconducting proton phase, or, conversely, determine the mass of Cas A from a reliable 
calculation of the size of the proton superconducting phase in superdense neutron star matter. We 
emphasize that the strategy outlined in this paper can be applied to any other neutron star of similar 
age, whose temperature might be reliably monitored over a several years period. High-mass neutron 
stars, such as the recently discovered neutron stars J1614-2230 (1.97 ± 0.04 Mq) and J0348-I-0432 
(2.01 ±0.04 Mq), appear particularly appealing as a significant fraction of the protons in their cores 
may be superconducting. 

PACS numbers: 04.40. Dg; 21.65. Cd; 26.60.-c; 97.60.Jd 



Introduction. One of the most prominent ways of 
probing the poorly known core compositions and inner 
structures of neutron stars consists in studying the ther- 
mal evolution of such objects theoretically jll-l23| and 
then comparing the outcome with the ever-increasing 
body of observed data on neutron star temperatures 
[1, [IS HI] . The quality of observed neutron star 
temperatures has culminated in the recent temperature 
data obtained for the compact object in the supernova 
remnant Cassiopeia A over a 10-year period |24{ . The 
observed temperature evolution of Cas A exhibited a 
steeper temperature drop than originally thought pos- 
sible for a star of that age (~ 300 years). Recent in- 
dependent studies have shown, however, that the ob- 
served behavior can be explained if the protons form a 
^5*0 superconductor with a pairing gap that is large and 
non- vanishing throughout the entire stellar core, and the 
neutrons in the core form a.^P2 su per fluid with a crit- 
ical temperature of 5 x 10^ K [I^, Uli- The relatively 
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sharp decline of the temperature of Cas A is then caused 
by the emission of neutrinos from neutron pair-breaking- 
formation processes [T^, [HI ■ It has also been shown that 
in-medium effects play a very important role for the tem- 
perature evolution of the neutron star in Cas A [2^ , and 
that the late onset of the direct Urea process, driven by 
the changing core compression during stellar spin-down, 
may provide an explanation of the observed data as well 

m 

Building on earlier work [lol [ll| , here we too assume 
that the neutrons in the core of the neutron star in Cas A 
have formed a '^P2 neutron superfluid. In contrast to the 
aforementioned work, however, we treat proton pairing in 
neutron stars [l^l on a more general footing by consid- 
ering a variety of different proton superconducting pro- 
files and critical temperatures. As shown in this letter, 
this allows us to establish an important constraint on the 
core density (core penetration depth) up to which pro- 
tons need to be superconducting so that agreement with 
the temperature data of Cas A is achieved. As in [lol.[Tl| 
we assume that the protons in the stellar core are paired 
in the ^Sq state since the first moments of the star's evo- 
lution, which implies a critical temperature of ^ 1Q'^~^'^K 
for proton superconductivity. The proton pairing is nec- 
essary as it suppresses both the (density dependent) di- 
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rect Urea process among nucleons, n p + e + Pe, as well 
as the modified Urea proeess, n + n n + p + e + v^, 
which would otherwise spoil the agreement with the Cas 
A data. 

Proton Superconductivity Depth. Wc assume that the 
matter in the core of a neutron star is well described 
by the Akmall-Pandharipande-Ravenhall (APR) model 
for the nuclear equation of state [2^, [2^ . We shall use 
a parametrized version of this equation of state, as de- 
scribed in [2^. The choice of parameters used in our 
study is the same as in [2^ , such that the direct Urea pro- 
eess is allowed for stars with masses greater than 1.4 Mq. 
Following [iH we allow the neutrons in the core of the 
neutron star to form pairs in spin-triplet states {'^P2)- As 
discussed in neutrons in the crust are likely to form 
spin-singlet (^S'o) pairs, but these have little effect on 
the thermal evolution. For the neutron triplet pairs in 
the core we use a phenomenological model, as described 
in p^ . The values of the supcrfluid gap, suppression fac- 
tors for neutrino emissivities, and the luminosity of the 
jjai r-breaking- formation processes were also taken from 

To determine how deep superconducting protons actu- 
ally penetrate into the core of the neutron star in Cas 
A, we start our cooling simulations for a superconduct- 
ing proton phase that penetrates deeply into the center 
of star. i.e. in this case all the protons are paired and 
have a very high critical temperature, meaning that the 
star begins its thermal evolution with all proto ns p aired 
up. As already known from the work of [13, [ill, this 
assumption leads to an agreement with the temperature 
data observed for Cas A. We then repeat the cooling 
simulation, but for incrementally reduced depths of the 
proton superconducting phase. This scheme is repeated 
until agreement with the Cas A temperatures is no longer 
achieved. The cooling simulations are carried out for neu- 
tron stars with masses of 1.4 M©, 1.6 Mq, 1.7 Mq, and 
1.8 Mq. As in [I^, all our cooling simulations are for an 
accreted neutron star envelope. The value chosen for the 
accreted mass is AM/M = 5.5 x 10"^^. 

In Fig. [1] we shown the temperature evolutions for 
our sample neutron stars, each one possessing a super- 
conducting proton core of proper size so that agreement 
with the Cas A data is achieved. The insert in the upper 
right corner is an enlargements of the temperature data 
of Cas A observe over a 10-year time period. Also show 
in Fig. [1] is the temperature evolution of a 1.8 A/q neu- 
tron star when all neutrons and protons in the core arc 
unpaired. The latter badly fails to describe the Cas A 
data. 

Figure [2] displays the critical temperatures of super- 
conducting protons for the neutron stars shown in Fig. 
[TJ Their central densities can be read off from the 
termination points of each curve. Figure [3] shows the 
density dependence of the ^Sq proton pairing gap for 
each stellar model. This figure shows that the depth 
of the superconducting proton phase increases with neu- 
tron star mass. This is a consequence of the direct Urea 
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FIG. 1: (Color online) Thermal evolution of neutron stars 
with different gravitational masses. The proton superconduc- 
tivity depth inside of each star has been chosen such that 
agreement with the Cas A data is achieved. (See text for 
more details.) 
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FIG. 2: (Color online) Critical temperatures, Tc, of supercon- 
ducting ^So protons as a function of baryon number density, 
n, for the stellar models of Fig. [1] (no — 0.16 fm^"^ denotes 
the baryon number density of ordinary nuclear matter.) 



process, whose fast cooling rate needs to be suppressed 
by the superconducting protons. Otherwise the stel- 
lar model would cool down much faster than observed 
for Cas A. Aside from the direct Urea process, super- 
conducting protons are also essential for controlling the 
changes in the stellar temperature caused by the modi- 
fied Urea process. It if were not for the superconducting 
protons, the core temperature would drop too quickly 
so that a confiict with the emission of neutrinos from 
supcrfluid pair-breaking neutrons would arise, again ren- 
dering agreement with the Cas A data impossible. In 
summary, we conclude that, in agreement with [lol. [ll|. 
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FIG. 3; (Color online) Gaps of superconducting ^Sq protons 
for the neutron star models shown in Fig. [21 Solid lines de- 
note a temperature of 10^" K, dashed lines correspond to a 
temperature of lO'' K. 
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FIG. 4: (Color online) Proton superconductivity depth as a 
function of neutron star mass. The solid line represents the 
best (2nd order) polynomial fit of the data. 

proton superconductivity needs to be present in the core 
of the neutron star in Cas A. Agreement with the ob- 
served temperatures of this object can however achieved 
for a range of different proton superconducting profiles 
and critical temperatures. The latter allows us to estab- 
lish a connection between the mass of a neutron star and 
the core penetration depth of superconducting protons. 



The result is shown in Fig. |4l where we plot the density 
(core penetration depth) up to which superconducting 
protons exist in the core of Cas A, considering a range of 
possible stellar mass. The solid line represents a 2nd or- 
der polynomial best fit to the numerical data. One sees 
that for lower stellar masses the curve flattens consid- 
erably, approaching a penetration depth of Pc/ po — 3.75 
for a canonical neutron star mass of 1.4 Mq. The flatten- 
ing is consistent with the fact that for low stellar masses 
the direct Urea process is highly suppressed due to low 
proton fractions inside of such stars. Hence, the super- 
conducting protons only need to suppress the modified 
Urea process along with other less rapid stellar cooling 
processes. This relationship shown in FigU] can be used 
to cither predict the depth of the superconducting proton 
phase, or, conversely, determine the mass of Cas A from 
a reUable calculation of the size of the proton supercon- 
ducting phase in superdense neutron star matter. 

Summary and Conclusions In this letter, we have thor- 
oughly investigated the role of proton superconductivity 
for the cooling of neutron stars, building on earlier work 
of Our results confirm that superconducting 

protons are required in order to explain the recent tem- 
perature data observed for Cas A. Furthermore our re- 
sults agree with the maximum critical temperature for 
supcrfluid neutron established in [lol. [Tl|. 

In contrast to earlier work, however, we study a se- 
ries of phenomenological proton-pairing models to de- 
termine how deep into the stellar core superconducting 
protons actually penetrate. This allows us to establish 
a heretofore unknown relationship between the mass of 
the neutron star in Cas A and the penetration depth of 
the superconducting proton phase. This relationship can 
be used to either predict the depth of the superconduct- 
ing proton phase, or, conversely, determine the mass of 
Cas A from a reliable calculation of the size of the pro- 
ton superconducting phase in superdense neutron star 
matter. We emphasize that the strategy outlined in this 
paper can be applied to any other neutron star of simi- 
lar age, whose temperature has been reliably monitored 
over a several years period. High-mass neutron stars, 
such as the recently discovered neutron stars J1614-2230 
0-97 ± 0.04 Mq) [Hi and J0348-h0432 (2.01 ± 0.04 A/q) 
[5J], appear particularly appealing as a significant frac- 
tion of the protons in their cores may be superconducting. 

This research clearly underlines the need for improved, 
state-of-the-art calculations focusing on the role of pro- 
tons (as well as neutrons) in superdense neutron star 
matter. 
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